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We derive the sum rule for the spectral function of the stress-energy tensor in the bulk (uniform 
dilatation) channel in a general class of strongly coupled field theories. This class includes theories 
holographically dual to a theory of gravity coupled to a single scalar field, representing the operator 
of the scale anomaly. In the limit when the operator becomes marginal, the sum rule coincides 
with that in QCD. Using the holographic model, we verify explicitly the cancellation between large 
and small frequency contributions to the spectral integral required to satisfy the sum rule in such 
QCD-like theories. 
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S , I. INTRODUCTION 

' The recent discovery of strongly coupled quark-gluon plasma at RHIC !lM4| has spurred a significant theoretical 
effort towards understanding the properties of matter described by strongly coupled quantum field theories at finite 
temperature. The transport properties of such theories are much more sensitive to the strength of the coupling than 
stationary thermodynamic properties. In particular, the near perfect fluidity in such theories is vievifed as a tell-tale 
sign of the strong coupling. 

' fH ' • The transport properties of a theory are closely related to the spectral functions of the operators such as stress-energy 
^ I tensor. In particular, viscosity can be determined from the low frequency limit of the spectral function by the well- 
' known Kubo formula. The first-principles calculation of the spectral functions using lattice Monte Carlo techniques is 
a challenging task 043; especially in the low- frequency regime, and the prior knowledge of the properties of spectral 
functions is essential. Therefore, constraints on the spectral functions in the form of the sum rules for the integral of 
I— the spectral function has been the subject of recent attention [l3|4lf 

The discovery of the AdS/CFT holographic corresDondence |16l4l8| h as opened new possibilities to study thermo- 
' dynamics and transport in strongly coupled theories (for reviews see 19l425l|). Although not yet directly applicable to 
QCD, at least until we know its holographic dual, these methods allow one to study generic properties of the strongly 
coupled plasmas using model theories, such as A/" = 4 SUSY YM, or holographic models which incorporate QCD-like 
0^ ' features such as confinement. In particular, the shear channel sum rule derived in [l5| has been verified in the J\f = 4 
O SUSY YM using AdS/CFT correspondence. 

In this paper, we concentrate on the sum rule for the spectral function in the bulk channel, corresponding to 
uniform dilatation or isotropic expansion. The bulk channel sum-rule is trivial in a conformal theory such as A/" = 4 
^-H ' SUSY YM. However, nontrivial sum rules in QCD do exists and have been a subject of recent studies (l3l - [T5| . Unlike 
, the shear channel, the bulk channel correlation function (and associated bulk viscosity) is sensitive not only to the 
' strength of the coupling, but also to the amount of the scaling violation. Therefore we consider a holographic model 
with the simplest mechanism of conformality violation, which is similar to QCD. We shall assume that this model 
describes field theories where scale invariance is broken by the presence of a scalar operator O with dimension A_|_ < 4 
^ ' in the action. The scale anomaly, T^, is proportional to this operator, and in QCD this role is played by the gluon 
- - condensate operator G^yG^" / as- 

Thermodynamics and transport in such theories have been first studied in Ref. [26l - l28| . In particular, it has been 
shown that the speed of sound Cs = \/ de/dp approaches the conformal value 1 / \/3 universally from below [29l . Isoj . 

The bulk sum rule has been also tested recently in Ref. in a theory whose holographic dual is a dilaton gravity 
with the Chamblin-Reall dilaton potential [sl], which has the virtue of being analytically tractable. However, the 
mechanism of the scaling violation in such a putative field-theory (if the field-theory dual exists) would differ from 
that in QCD. In particular, the speed of sound does not approach its conformal value in the high-temperature limit. 

In the general class of theories we consider, we shall find that the temperature dependent part of the bulk spectral 
function [^(w)]^, obeys the following sum rule 

34-A.)(.-3rt = irM^<^. ,1) 



ds I TT ./n W 
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In the limit when the operator O becomes marginal, A+ — 4, the sum rule coincides with the sum rule in QCD, 
derived in [l^ . 

We shall also find that in this lirnit, which we shall refer to as marginal, the sum rule in Eq. ([T|) exhibits the same 
puzzle as discussed in Refs.fl^ [s^ [sj in the context of QCD: the l.h.s. of the sum rule is of order a^, while p ~ a^, 
where is the QCD coupling. Here, in the marginal limit, the l.h.s. of Eq.(l) is of order A'i while p ^ A?_, where 

A_=4-A+. (2) 

We shall demonstrate that a delicate cancellation indeed occurs between the high frequency tail of the integral, 
uj ^ T, and the intermediate region of w ^ T, of the kind needed to occur in QCD. Keeping only the leading 0{A'^) 
terms in [p\j., we find that the integral in Eq. ([l} converges. However, the integral over the contribution of the 
subleading 0{A^) terms in [p]j, has support over an interval of w which stretches to infinity as 1/A_ in the marginal 
limit A_ 0. We evaluate the resulting additional 0(A?_) contribution from this long tail to the r.h.s. of Eq. ([T]) 
analytically and show that it cancels the O(A^) contribution from the lo ^ T region. 

This paper is organized as follows. In the next section, we present the definitions of quantities involved in the sum 
rule and its derivation. In Section Hill we derive the sum rule, Eq. ([1]), in the general class of theories we consider. 
Section ITVl introduces the holographic description of such theories. Section Ivl explains our method of calculating the 
spectral function and the related Green's functions. In Section I VII we make analytical calculation of the large 
asymptotics of p{uj), which we use in Section IVHl to subtract from the numerically determined p{uj) to obtain 
The sum rule is verified numerically in Section rVIIII for a sample of values of A_. We demonstrate the cancellation 
required to satisfy the sum rule in the marginal limit in Section IVIIII Cross-check of the results with the existing 
analytical result for bulk viscosity (Ref.fss]) is made in Section [Xj We conclude in Section IxH Appendix \K\ contains 
the relevant results from [29,] , used throughout the paper. 



II. DEFINITIONS 

We shall consider only the response to homogeneous (q = 0) perturbations and define the spectral function for the 
trace of the stress-energy tensor T'"^, as usual, by 

p[ijj) = — ImGfl,(a;), with (3) 
Gr{^) ^ -I / die-* / d^x {m{x),T:m)- (4) 



The definition of the retarded Green's function Gr is subject to the usual ambiguity due to the contribution of the 
product of the operators at the same point at x = 0. However, these contact terms do not have imaginary parts and 
do not contribute to p. 

It is convenient, as in to define Green's functions of T^^^ by considering the response of the system to the 
perturbation of the background metric g^v = rj^i, + 5g^i, around the fiat Minkowski metric rj^^, = diag(— 1, 1, 1, 1). 
Euclidean time correlators can be obtained similarly from variations of the Euclidean partition function Z[g]: 

\ogZ[g] = \ j d^x{T^^'^{x))6g^,{x)+^- j j d^xd^y{T^^'^{x)T^P{y))59,,,{x)5g,^p{y^ (5) 

The correlation functions of the trace 9 = can be also defined via variations of the partition function with respect 
to the metric variations of a special form (dilatations): 

gmt.. - '7M.e-2f^ (6) 
\ogZ[gm = - fd^x {e{x)Mx) + 1 f f d^xd^y {9{x)9{y)) n{x)n{y) + (7) 



2 _ _ 

The two definitions of the two-point function of the trace differ by a contact term: 

{9{x)) = v^AT^'i^)) (8) 
{9{x)9{y)) = rj^,r,xp{T^'^{x)T>^''{y)) + 25\x - 2/)r?^.(T'-(x)) (9) 

The corresponding retarded correlators can be defined via linear response to perturbation VL: 

d 



Gr{x, y) = {9{x)9{y))R = -9{^{y)my))) (10) 
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The definition (ITOl) is convenient because the value of its Fourier transform (for q = 0) 

Gr{u:) = {00)r{lu) = fdte''^' f d^x {e{x)e{0)) r (11) 



at vanishing frequency oj fohows from the conservation of entropy in the ideal hydrodynamics [l5j | : 

mRiO) ^ ^ [,/^){e)) ^ (^3s^ - 4) {6} (12) 



where s is the entropy density. 



III. THE SUM RULE 



The sum rule derived in jlSj applies to the zero temperature subtracted Green's function 

[GrIt = Gr, ~ g'^j^ (13) 
and its imaginary part on the real axis [p\rp — — Im [Gr]^^: 

[Gr{ioo)]t - [GrMt = - Hdu lPMk_ (14) 

Using Eq. (IT^ and [{0)]rp = 3p — e, one can write 

[Gi?(0)]^ = - (^3s^-4^ (6-3p) (15) 

where e = [(T™)]^ and p = [{T^^)]j, are equilibrium thermal energy and pressure at given entropy density s. In 
QCD, due to the asymptotic freedom, [G R{ioo)]rp — and one obtains the sum rule found in fisj : 



33— ~ 4 (e - 3p) ^ - \di^ ^'^ (16) 
OS J n J LO 

In this paper, we consider a generic conformal field theory perturbed by an operator O of dimension A+ sourced 
by the field c. We place this theory on a nontrivial gravitational background given by metric g^^. The change of the 
partition function of the theory under dilatations © is equivalent to the rescaling of the only dimensionful external 
field c, the source of the operator O: 

Z[g(f7),c] = Z[5(0),e-^-"c], (17) 

where A_ = 4 - A+. Thus 
Thus 



A_e-^-''c^^ = A_e-^-^'c(0). (18) 
oc 

{e{x))^-A.c{0{x)), (19) 



and 

{eix)9iy)) = A^_c^{Oix)0{y)) - A^_c{0)5\x ~ y). (20) 

To evaluate the Fourier transform of the correlation function (j20p at large (imaginary) frequency w, we can use 
the operator product expansion (OPE) for 0{x)0{y). The leading contribution comes from the unit operator, and 
behaves as but cancels in [Gi^]y because it is independent of temperature. The contribution of an operator 

of dimension A comes with the Wilson coefficient which behaves as uj'^^+~'^~^ . Assuming that there are no operators 
(with vacuum quantum numbers) of dimension equal to or lower than 2A4- — 4, we conclude that the contribution of 
the first term in Eq. (|20l) vanishes in the u — )■ ioo limit. The assumption implies, in particular, that A+ < 4, i.e., 
A_ > 0. The contribution of the last term in Eq. (PH)) then gives the value for GR{ioo): 

[Gfl(z(X))]j, = -[{e0){i^)]^ = A^ c [{0)]j. = A_ (e - 3p), (21) 

where we used the fact that analytical continuation of a retarded correlation function to Matsubara frequencies on 
the imaginary axis equals negative of the Euclidean correlator. We shall verify Eq. (1211) explicitly and analytically in 
the holographic model. Combining Eqs. (|15p . 1^1^ and we obtain Eq. ([T|). 
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IV. HOLOGRAPHIC MODEL 



As discussed above, we will consider a four-dimensional (4D) conformal theory in which the conformal symmetry 
is broken by the operator O with scaling dimension A+. The holographically dual description of such a theory 
must therefore include a five-dimensional (5D) scalar field with mass m| = — A_|.A_. As a model for the QCD 
thermodynamics, such a theory was first considered in [26l [27| . It models in the most straightforward way the 
breaking of the scaling invariance in QCD. The minimum set of fields required to describe correlators of the stress 
energy tensor and of the operator O include the 5D metric gMN {M, N = 0,1, 2, 3, z) and the scalar field (dilaton) (j). 
The minimal action is thus given by 





(22) 



where R is the Ricci scalar, g is the determinant of the metric, 7 is the determinant of the induced metric on the UV 
boundary dM a,t z = 0, K is the extrinsic curvature on dM, and is the 5D Einstein gravitational constant. The 
value of is inversely proportional to the number of the degrees of freedom in the dual four-dimensional theory, e.g., 
in a gauge theory with large number of colors Nc- The smallness of (i.e., the largeness of the number of colors) 
controls the semiclassical approximation which we use. The last term in the action is the Gibbons-Hawking term, 
which removes the boundary terms arising upon integration by parts of the terms in R linear in second derivatives of 
the metric [36| . 

The potential for the dilaton, V{(l)), is the function which was tuned in [2^[23| to best "mimic" the QCD equation of 
state. Here, similarly to |2^, we shall concentrate on the results which are universal in the class of models described by 
Eq. (22) with any (sensible) potential. For example, Ref. [2^ found that the speed of sound approaches the conformal 
value l/VS as T — >■ 00 universally from below in such theories. Here we shall also make use of the large T limit to 
the extent that it makes only the curvature of the potential V"{(p) = m\ = A+ (A+ — 4) matter. 

The negative cosmological constant provided by choosing V{0) = — 12 ensures necessary asymptotics of the metric 
near the boundary gMN ^ while the dilaton field asymptotics is fixed by 1^"(0): (f> ~ z'^- . The singular behavior 
at z = can be regulated by setting the boundary conditions at z = e and taking the UV regulator e to after 
necessary renormalization. The rules of the holographic correspondence require us to extremize the action S'5 w.r.t. 
the metric and the dilaton field subject to the boundary conditions on the UV boundary z — e —?' 0: 

9iiy{x,z)\^^^^ g^,^{x)£^'^, (j){x,z)\^^^ = c{x)£^- . (23) 

The holographic duality then implies that the extremal value of the Euclidean 5D action as a functional of the 
boundary values g^u{x) and c{x) is the same functional as — log Z[g^u,c\ in the 4D theory. In order to calculate the 
Green's function {69) we use Eq. © and the relationship ^ between the derivatives of log Z and the correlators of 

The one-point functions {T^'^) are determined by the solutions of the equations of motion with boundary conditions 
homogeneous in the 4D coordinates x''^. We follow Ref. |2^ to determine the solution to the equations of motion 
homogeneous in 4D. In the gauge chosen in Ref. (29j . the metric has the form 

ds^ = ^ (^-/ (z) dt^ + dx' + e'^'^'^J^) ^ g^dx^dx"", (24) 

where f{z) and B(z) are functions only of the extra dimensional coordinate z, and f{z) has a simple zero at some 
value of z = z//. The functions, / and B, and the background dilaton field, 0, are all determined by extremizing the 
action. The details can be found in and for completeness presented in Appendix [Xj 

The extremum of the action with boundary conditions (j23p is a one-parameter family of solutions to the equations of 
motion. A convenient choice of the parameter is the enthalpy density, w, which arises as an integration constant (see 
Eq. (jA6[) ). The temperature is determined, as usual, by considering the periodicity in the Euclidean time necessary 
to avoid conical singularity ai z = zh- 

r=^e-^(^")|/'(z^)|, (25) 



and is related to w via equation of state (see e.g. Eq. (jAlSp ). 
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V. TWO-POINT FUNCTIONS 



The calculation of two-point correlation functions requires us to consider non-homogeneous solutions to the equa- 
tions of motion. We only need to consider infinitesimal variations around the homogeneous solution (j24p and expand 
the action to quadratic order in these variations. We parameterize these inhomogeneous solutions as 

9m N = 5MAr(l + It-mn) (26) 

where summation over AI, N is not implied. We exercise the freedom of the gauge choice to set z) = (/)^^\z), i.e., 
we set variations of around homogeneous solution to (we can do this as long as we are not interested in calculating 
correlation functions of O, i.e., as long as c remains homogeneous). We find that this gauge choice provides nontrivial 
simplification of the equations of motion, allowing us to reduce them to a second order equation for metric variation 
H, instead of the third. 

We then consider solution homogeneous in the spatial coordinates a;*, and depending only on a;*^ = t and z, since we 
are interested in q = variations of 5^^. For such solutions we can use the remaining gauge freedom to set hzi_i = 0, 
but not hzz- 

We can also use 0(3) symmetry to simplify our analysis by separating the spatial part of metric perturbations into 
trace and traceless parts: 

h,j^H5,^+hl. (27) 

At quadratic order, H mixes with components Hqq and hzz, but decouples from the traceless part hj^ as well as from 
the off-diagonal components /loi- Therefore, we shall focus only on the metric perturbations H, hoo and hzz- We can 
express the coupling of the metric perturbation to T'^'^ as 

h^.T'^'' = i/S + hf^T^' + 2h^,T°' + /looT"", (28) 

where E = T'^^Sij is the trace and = T'-* — is the traceless part of the stress tensor T'^. The two-point 

function {96) in Eq. ([9]), since r/i^i^T'^^ = E — r"°, can be expressed in terms of correlation functions of and E 
using 

V,..VXp{T^''{x)T''iy)) = (T°°(x)T°"(y)) - 2(r°"(x)E(y)) + (E(.T)E(y)). (29) 
The corresponding retarded correlators can be found using holographic correspondence and the recipe (37l |: 

5^ ^5 



(T""'(t)T""'(i')) 



R 



Shooit, z)ShoQ{t' , z) 



(T™(t)E(<'))7?, = -4 



5hoo(t,z)5H{t',z) 
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SH{t,z)SH{t',z) 



(30) 
(31) 
(32) 



Since we need the Fourier transform of (99) it would be convenient to express 5*5 directly in terms of the Fourier 



hMNit,z) = I ^hMNiuJ,z)e'"^K (33) 

ZTT 



modes of h^^i, defined as 



The metric perturbations H, hzz and hoo, defined by Eqs. (|26p and (|27|) . satisfy a set of equations derived from the 
linearized Einstein's equations: 



hzzioJ,z) = -zH'{lu,z)+H{uj,z)(^^ 



z f{z) 
fiz) 



(36) 
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with the primes denoting derivatives with respect to z. In order for Eqs. ([30|) - (|32)) to give the retarded correlation 
functions, the solutions to Eqs. must satisfy the in- falling wave conditions at the horizon z = zh- 

In order to calculate the two-point functions, we expand the holographic action in Eq. ([22l) to quadratic order in 
metric variations around the homogeneous background. This expansion can be written in a compact form; the bulk 
part of the action is given by. 



'5'bu; 



1 f duj 



Ik 



2k2 / 27r 



■ d X dz ■ 



'B(z) 



'H"HhH + H'Hl2H' + H'^fMsH + HHd^H + (c.c.)) , (37) 



where the dagger refers to the transposed complex conjugate {H*{uj, z) — H{—uj, z)), and 



H 



H{io,z) 
hoo{uj,z) 
hzz{^, z) 



3 -3 -3 
-3 -1 1 
-3 1 3 



-3 3 3 
Ml = ^ I 3 1 -1 


-, / 3 -3 -3 

-B'{z) I -3 -1 1 




3 3 
Ma = - I 3 2 -1 

3-10 



4 f{z) 



-6 6 6 
9 3-3 
3 -1 -3 



3 12 \ 3^2g2B(.) / 2 1 

12 -4 -9 / ^ Jy^^) \ 1 



(38) 



Integrating by parts and using equation of motion, we reduce the action at the extremum to the boundary term: 



ShuW 



1 f doj 



2k^ / 2tt 



■d^x- 



-B{z] 



fiz) 



\f{z) z 



2H'^M2H' + HH'hH 



(39) 



The action also receives a contribution from the Gibbons-Hawking boundary term: 



1 fduj ^ 



'B{z) 



fiz) 



{h'^^M^H + H^'MeH + (c.c. 



with 



Ms = Ah and Ah = ^ f- - 

8 f[z) 



3 -3 -3 
-3 -1 1 
-3 1 3 



(40) 



Combining all of this, the action at the extremum can be expressed in terms of the boundary values of hoo, H and 
H': 



S5 



1 f duj 



-B{z) 



2k2 / 2tt 



T^d-'x- 



fiz) 



—K„{uj, z)hoo{uj, z) + -— -I- - 



2z 



9 , 3/'(z) 



2z 4 f{z) 



{Ka{uj,z)H{u,z) + {c.c.)) 



3 B'{z)f'{z) 3 f'{z)^ 3 uj^e^^<-'^ 
2^2 + 4 f{z) 8/(z)2 2 /(z)2 



zH*{uj,z)H{uj,z) 
- ^zB'{z)H*{io,z)H'iL,,z)^ 



(41) 



From this one can immediately read off the expressions for the Fourier transforms of the stress-energy two-point 
functions: 



(r"°T°")(w) - 



1 6e-'SW/(z) 



2k2 



00 \ 



-{T 



6e-«(^)/(z) e--»(^)/'(z) 



2 B'{z)f{z) , f(z)2 2a;2g2i?(^) ^ 2B'{z) H'jio, z) 



3/(z) 6/(z)2 Sfizy 



3 i/(w,z) 



(42) 
(43) 
(44) 
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where in the last equality in Eqs. (|42)) and (|43|) we used expressions for one-point functions Ref. {29j. As ex- 
pected, the correlators involving T^^ are frequency independent and satisfy energy-momentum conservation Ward 
identies Ref. [sMi- The above expressions contain contributions divergent in the limit e ^ 0. They can be removed 
by subtracting a polynomial of sufficient to cancel the divergences. These terms do not affect p{io). 
Putting this all together, the finite (as e — >■ 0) part of the bulk Green's function is found to be. 



GRiiu) - P{u;^) = 4(e - 3p) + 



6 e-'^(^)/(^)S'(z) i/'(cj,z) 



2k2 



H{uj,z) 



(45) 



where the polynomial P{ijJ^) = 4(T™ — I])t=o+w^ •Se^*-^-' / (re^e^) combines the UV divergent, temperature independent 
contact terms. The bulk spectral function is then given by 



p{uj) 



6 e-^'-^\f{z)B'{z) H'{uj,z) 
2^ ^2 H{uj,z) 



(46) 



VI. LARGE FREQUENCY ASYMPTOTICS 



In order to calculate the bulk retarded Green's function and the bulk spectral function, one needs to solve the 
equation of motion for the field H{uj,z), Eq. In general, this cannot be done analytically. However, in the large 

w limit, the solution can be found similarly to the Born approximation in quantum mechanics. 

We perform a Louiville transformation to bring Eq. p4p to the Schrodinger form. In terms of the new coordinate 
X, given by 



and the new function 'I'(x), given by 



Eq. ([M]) takes the form 



7(^^"'' 
H{z) = zB'{z)-^/H{x), 



with the Schrodinger potential, as an implicit function of x, is given by 



Uschix) 



f{zf / 2 1 B"{z) 1 B"{zf 1 B"'{z) 1 f'{z) 
e'^Biz) [ji ^ ~ B'{z) ^ 4 B'{zY ^ 2 B'{z) ^ Yz f{z) 



+ Vz)-V(z) 
fiz) B'{z)+ z""^'^ 2^' 



(47) 



(48) 



(49) 



(50) 



The in-falling boundary condition on H{z) at the horizon z — zh transforms into the outgoing wave condition on 
^(x) at X = c». The bulk spectral function can be expressed in terms of the new function 5' as 



P 



6B'(z) ^ ^'(x) 
t; — Im- 



2k2 z2 ^/{x) 

or somewhat more intuitively (using Eq. (jAip ) as 

0'(z)2 Im[**(x)*'(x)] 



2k2z 



|*(X)|2 



(51) 



(52) 



which has the quantum mechanical interpretation of the probability flux for the wave function (normalized as 
|5'(x)|2 — 0'(e)2/(2K2£) at z = e). The positivity of the spectral function can be seen then as a consequence of the 
conservation of the flux and the outgoing wave boundary condition at x = oo. 

At large w, we would like to treat the potential as a perturbation and calculate the "wave function" 'J using the 
Born approximation. There is one difficulty, however: the Schrodinger potential Ugchix) diverges at a; = 0. We shall 
separate the leading divergence explicitly 



(f^ix) ,,,(2 (3-2A_)(2A+-3) „^ 

+ ^(x) a;2 - ^ '\ '- - SUschix) 



(ix2 



4x2 



= 



(53) 
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and treat the remaining part of the potential, 5Usch{x), as a perturbation. We can calculate ^'(x) iteratively, using the 
Green's function method. The solution to Eq. (1531) with SUsch — 0, satisfying the outgoing wave boundary condition 
at a; = oo, is given by, up to unimportant normalization, 

^oix)^iLuxy/^Hl^\uJx) (54) 

where v — 2 ~ A_, iJ^^-* is the Hankel function the first kind. Substituting ^I^o into P5|) and then into (H5|) we find, 
using Eq. (|X9)) for B'{z), 

a„M - P(^') - ^f^l^ (»t(.A_) - ,:, {^f^-"- -^Al (A, - A_) + (55) 

where the dots stand for corrections to the leading order which would come from iterations of 6Usch- The first term 
on the r.h.s. is the leading term in the w — > cx) limit. It is temperature independent and has correct w-scaling to be 
identified with the leading contribution of the unit operator to the OPE of (OO) in Eq. (pn)) . This term has nontrivial, 
but temperature independent, imaginary part and it is subtracted when we calculate [G^jJ-^. 
The second term on the r.h.s. of Eq. (|55|) gives 

[Gfl(*oo)]j, = (A+ - A_) = A_(e - 3p). (56) 

where we used Eq. (jAll[) . This agrees with Eq. ([2T|) as expected. 

The leading correction, ^ti, to ^Pq can be found, similarly to the Born approximation, by using Green's function 
G{x,x'): 



^lix) = / dx'Gix,x')^oix')SUsch{x'). (57) 



The Green's function satisfies equation 



d^Gix,x') ^. / 2 (2A_ - 3) (2A+ - 3) 



+G{x,x')[c.'- ' - " '-j=s{x~x'), (58) 

with boundary condition G{x, x') — ai x — e and outgoing wave condition at a; = oo. The solution is given by 

, J -/<ix') (a;a;)i/2j^(tjx) x < x' 



where 



l<ix')^~—{u:x'y^'Hl'\u:x') and j^{x') = -— {lox')'/' Mojx') . (60) 
2uj 2uj 



Substituting into Eq. ([57)) wc find 



x'SUschix') (hI^^ujx')) dx'+Hl^^ujx) / x' 6Uschix')Hl^Hujx')J^iu}x')dx' 



*i(x) = -^(c..)V2 



To zeroth order in 6Usch, \I/ = vE'o and Eq. (|5T|) gives 



(61) 



2K2r(2- A_)2 V2 
To next order in SUsch, ^' = 5*0 + and p ~ po + Pi, where 



PoM = TTIoT^T^^-T^ o ■ (62) 



pi(w) — 7rpo(w) / X SUsch{x)J^{(jJx)Y^{ujx)dx . (63) 



In the limit a; — >■ oo, the integral in Eq. ()63p is dominated by the region of small x ^ u! ^. Expanding SUsch in 
powers of x, we generate a l/w series (possibly asymptotic) for the integral. 
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To enable further analytic calculation, we shall limit it to leading terms in the expansion in powers of (? . As we 
also discuss in Appendix 1X1 this corresponds to the regime of high temperatures, i.e., T c^/^-, since c is the 
only dimensionful parameter in the theory. In addition to enabling analytic calculations, this limit has the virtue 
of yielding results which do not depend on the form of the dilaton potential beyond its curvature at the minimum 
(related to A_). In this regime, to the order we work, we can set c to zero in J7sch- The leading term in the Taylor 
expansion of t/sch is given by 

^U^.Y.ix) = ^ A_ A+u; + (64) 

where w = 2k^w, as defined in Appendix \F\ and 

Pi = -^^A_A+(1 - A_)(A+ - 1)(A+ - A_)po(^). (65) 

Combining Eq. ([62]) and Eq. ([65]) we see that p\ ~ w"^^- — >■ as w — >■ cxd. 

We note that, strictly speaking, not all O(c^) corrections to this result are negligible, since some of them grow with 
cj. However, as one can check, those corrections do not depend on w (i.e., on T) and are part of the ©(c^w^^'^- ) 
corrections to the Wilson coefficient of the unit operator in the OPE of (OO) . These terms would be also subtracted 
if \p(iS)\rj^ was calculated to the next order in (? . 



VII. CALCULATING THE SPECTRAL FUNCTION 



In order to calculate the bulk spectral function, one needs to first solve the equation of motion for the functions 
/(z), B{z) and (/)(z) and then use those results to solve Eq. (p4l) for H{z). This cannot be done analytically for an 
arbitrary potential V{4>). However, the first step of this calculation can be done analytically in the limit where the 
temperature is large compared with the conformal breaking scale, T ^ 0^/"^- . Another advantage of this limit is that 
the results are universal in the sense that the only property of the potential V{(p) which matters is the curvature at 
the minimum, which determines the dimension of the operator A+ . 

Equation would still need to be solved numerically, however. We use a shooting procedure starting from the 
IR boundary, i.e., the horizon z = zh- The equation has a regular singular point at z — zh ■ A Frobenius expansion 
around this point can be constructed iteratively. Among the two linearly independent solutions the in-falling wave is 
chosen. The expansion is used to calculate H and H' at a point near the IR boundary. These values are then used 
as the initial conditions to integrate the differential equation numerically all the way to the UV boundary z = e. The 
spectral function can then be determined from the numerical results using Eq. ()46|) . 

As a representative example, in Fig. l(a)[ we plot the spectral function at order for A_ — > 0. As expected, it 
is non-negative and it diverges for large w. The subtracted spec tral function [p{oj)]rp to order 0{c^) is obtained by 
subtracting po determined analytically in Eq. (|62|). Figures 1(b) and [2] show [p{u})]j,/uj plotted at different values of 
A_, while the dashed line corresponds to the analytic expressions for the large lu behavior given by pi in Eq. 

It is interesting to compare these plots with the lattice calculations of the spectral function in Fig. 5 of Ref. 
In both cases one finds the sign oscillation of the spectral function in the region lo ^ 2t:T, which is absent in the weak 
coupling (Boltzmann equation) result [H, [sl] . 



VIII. VERIFYING THE SUM RULE 



In the high-temperature limit T ^ c^^^- , the left-hand side (thermodynamic side) of the sum rule Eq. ([T]) can be 
analytically calculated in the holographic model. The trace anomaly is given by Eq. (jA14[) . while the derivative with 
respect to entropy can be expressed via derivative with respect to enthalpy w. Using the expression for d in the high 
T limit, Eq. (|A19|) . the derivative can be performed, and the left-hand side can be expressed as, 

{^'^s - ~ = ^ ^- ^-^ ■ (^^^ 

The right-hand side can be numerically calculated from [p]rp. For efficiency, we make use of the analytical result for 
the asymptotics of [p]rp Eq. We calculate numerically only the contribution to the integral up to some relatively 

large frequency a;,nax- The contribution from Wmax to oo is calculated analytically using asymptotics [p\rp w pi from 
Eq. ([65|) . We chose Wmax = IOttT. Table I] shows that the sum rule holds to at least within a fraction of a percent. 
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p(w)lu 




(a) (b) 



FIG. 1: The bulk spectral function divided by frequency in units of ^^c^ Ai {i^T)^ for A_ — >■ 0. a) The un-subtracted function, 
b) The function after the T = subtraction. 




FIG. 2: The bulk spectral function divided by frequency expressed in units of ^^c^ (ttT)^"^^" for A_ equaling 0.2, 0.5, 
0.7, and 1.0. The dashed curve in each plot is from the analytic expression for the spectral function in the large uj limit in 
Eq. (|65|) . The dotted line is drawn at the value of the bulk viscosity calculated from Eq. (I7ip in units of rj^^c^ Al (ttT)^^^^" . 
This value should be compared with the intercept of the plots at = 0. 



IX. THE MARGINAL LIMIT AND THE TAIL OF THE SPECTRAL FUNCTION 

As A- — > the left-hand side of the sum rule given by Eq. (|66l) scales as A'i [d ^ A_ according to Eq. (jA19P ). 
This is also apparent in Table H] 

However, the spectral function [p\rp on the right-hand side scales as IStl. This comes from factor B'{z) in Eq. p6)) 
which scales as A^ according to Eq. (|A9p . Therefore, in order to determine the ©(A^) part of p{ijj) we can simply 
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A_ 


LHS 


RHS 


% Error 


1.0 


-0.4569 


-0.4588 


0.4 


0.7 


-0.4862 


-0.4889 


0.6 


0.5 


-0.4951 


-0.4984 


0.6 


0.2 


-0.4997 


-0.5025 


0.5 



TABLE I: Verification of the sum rule in Eq. ^ in the holographic model for a sample of values of A_ in the high-temperature 
limit T » c^^'^- . The columns LHS and the RHS are the values of the left- and the right- hand side of Eq. ([l| expressed in 
units of 2^c^A?.(7rT)3-2'^-, which removes the dependence on k, c, T and the leading (A?.) dependence on A_. The LHS is 
calculated analytically using Eq. H66p . while the RHS is a numerical calculation of the integral as described in text. The last 
column shows the percent discrepancy between the two sides. 



set A_ = in Eq. ([34|) which then becomes 



H"i^)^H'i^)[l-^)--^Hi^). (67) 

where / is given by Eq. (jA16p in the high-temperature limit. This is the well known equation of motion of a massless 
scalar field on a pure AdS background. However, the analytic solution to this equation with required boundary 
conditions has not yet been found (see [iol - lisj for discussion). Nevertheless, the solution can be obtained numerically 
as d iscuss ed in Sec. I Villi The resulting spectral function is displayed in Fig.[TJ The O(A-i) part of [p{uj)]rp shown in 



Fig. 1(b) oscillates around 0, but the contribution from it to the spectral integral in the sum rule is non- vanishing. 



We calculate the integral numerically and find 

2 r-- [pHIt ^ 2c2a2_ f -) (.60000(1)) + 0(A3_) , (68) 

where w is the enthalpy. The integral of the A^ part of [p(a;)]/jn is converging very fast (apparently exponentially) at 
a; = oo and we picked relatively large Wmax = IOttT for this numerical calculation. The number in parenthesis is 3/5 
to at least four significant digits. How then is the sum rule satisfied if the left-hand side is only O(Ail)? 

We shall find the missing contribution in the tail of the function [^(a;)]^^. Indeed, taking Wmax to oo requires extra 
care. Because ©(A^) part of [p{uj)]j, decreases very fast as w — ^ oo, for sufficiently large u! the dominant part in 
[p(a;)]y is of order A'i. As we see analytically in Eqs. ((65)) . (l62t . [p]rp — >■ pi ~ A^uj~^^- . At any finite A_ the 
negative power-law tail effectively cuts off the integral at large cjtaii ^ exp(l/(2A_)). However, the contribution of 
this long tail to the integral grows with Wtaii as doj [p]rp/uj ^ log(a;taii/winax) ^ i-c, the length of the 

tail compensates for the extra power of A_. 

Let us calculate this 0(A?_) contribution from the tail. Using Eqs. (p5)) . ([5^ we can write for large uj and small A_: 

[p{uj)]^ =p,{uj) + ... = -^TTc'A'^u:-'^- ( j) + . . . . (69) 

Thus 

' '^[p],^-2c^A^(^)(l)+OiAl). (70) 



We see that the contribution of the high-frequency tail exactly cancels the contribution from the region oi u ^ T in 
Eq. ([55)1 . Therefore the RHS of the sum rule is proportional to A'^ just as the LHS. 

Remarkably, this cancellation mechanism is very similar to the one which was found in QCD by Caron-Huot in 
Ref ■ [34I . In the case of QCD the left-hand (thermodynamic) side of the sum rule is of order a^, while the spectral 
function on the right-hand side is 0{al), with logarithmically long high-frequency tail of order 0{al). 

We would like also to comment that the A?_ contribution to the bulk spectral function [p]j, in the high-temperature 
regime is the same, up to a constant, as the spectral function in the 4 SUSY YM theory. In fact, the spectral 
integral in Eq. (|68p has been also performed numerically in Ref. using a different method and with the same 



Note that according to Eq. I)A21|I . B'(z) ~ c^, which means, to O(c^) we are working, we can drop terms like B'{z), but not terms like 
B"(z)/B'{z) -1/2. 
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result. Unlike the case we consider, however, for the shear channel sum rule considered in Ref. [15| that integral 
saturated the sum rule. 

Also worth noting is that the shear and the bulk channel spectral functions are proportional to each other in the 
Chamblin-Real dilaton model examined in [3T| and "ii"] . In such theories the mechanism of the saturation of the bulk 
sum rule is the same as that of the shear sum rule and no high/low frequency cancellations are needed in either case. 



X. VISCOSITY 



As another cross-check of our results, we can compare them to the analytic calculation of the bulk viscosity in the 
high temperature limit of Ref. [35] , which finds 

1 C^Al . ^^3-2A_ .A_^ { r(l-A_/4) \ ' 

^=9^^"^^ 2 " U(l/2-A_/4)j • 

This formula was derived by matching gradient expansion of the stress-energy tensor to the gradient expansion of the 
background metric. 

Ideally, one should also be able to derive Eq. ([7T|) by solving Eq. at small w and using Kubo formula 

C - Ip'iO). (72) 



Indeed, this can be done analytically for A_ — 0. In this case, Eq. (|34l) reduces to Eq. ((67|) . For small w, this equation 
can be easily integrated. Normalizing as H{0) = 1 and using in-falling boundary condition at the horizon we find ^ 

Hiz) = 1 - !^ (I) log(l - wzV4) + .... (73) 
This gives for the spectral function at small w, 

p=^{7rTf^ + ..., (74) 

and for the bulk viscosity 

1 c^A^ 

('l^i-n'. (75) 

This agrees with the A_ ^ limit of Eq. ((7T|) . 

Unfortunately, a similar approach at finite A_ appears intractable, since an analytic solution to Eq. p4|) . even at 
small Lu, for arbitrary A_ is not known. However, the bulk viscosity can be calculated by applying the Kubo formula 
Eq. (17^ to our numerical results for the spectral function. We have verified the agreement of Eq. ([7T|) with these 
numerical calculations and illustrated it in Fig. [2] 



XI. CONCLUSION 



We studied the spectral function corresponding to time-dependent bulk deformation (uniform expansion) in a class 
of field theories where conformality is broken "softly" in the sense that at high temperature the equation of state 
approaches conformal limit e = 3p, similar to QCD. The breaking is due to a scalar operator of conformal dimension 
A+ < 4, which is an analogue of the operator of gluon condensate in QCD. We find that in such theories the bulk 
spectral function satisfies the sum rule given by Eq. ([!]). 

The sum rule in Eq. ([T|) is similar to the sum rule in an asymptotically free theory such as QCD derived by 
Romatschke and Son in Ref. . In the marginal limit A+ — > 4, the two sum rules are identical. We used this similarity 
to address an interesting puzzle noted in Ref. [33|]. In order to satisfy the sum rule, a delicate cancellation must occur 



^ Although this solution is not valid all the way to the horizon due to the singularity, for sufficiently small uj it is valid close enough to 
the horizon to allow matching it to the in-falling wave solution. 
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between the regions of w ~ T and of w ^ T in the spectral integral. Ref. [34| showed that this cancellation can be 
indeed seen in the high-temperature limit of QCD, using weak coupling calculation. Our study of the holographic 
model shows that such a cancellation is very generic to the whole class of strongly coupled theories with softly broken 
conformal symmetry. 

To make the connection of our model with QCD more tangible, we can observe that in QCD the action contains 
operator G^^yG^'^ /as and thus identify the operator O with G ^^G^^^ / as, up to a numerical coefficient. This coefficient 
can be determined by matching correlation functions (e.g., the bulk spectral function) of the holographic model to 
QCD, but we shall not need it here. The anomalous dimension of this operator in QCD is given by /3(as)/as = 
6oas/(27r). In QCD this anomalous dimension is a function of scale, vanishing logarithmically, i.e., very slowly, with 
increasing energy-momentum scale. In the class of QCD-like theories we consider, the corresponding quantity is 
A_|_ — 4 = — A_, which is a constant.^ This scaling dimension taken dX uj — oo, determines the value of [Gfi(ioo)]rp 
in the sum rule (|14l) , according to Eq. (I?!]) , and it is responsible for the difference of the general sum rule ([T|) from 
QCD. 

The holographic theories we consider can serve only as a qualitative or semi-quantitative guide to the QCD phe- 
nomena, since QCD becomes a weakly coupled theory at sufficiently high energy-momentum scale due to asymptotic 
freedom. But this guide might be useful by offering a view complementary to the weak coupling extrapolation out of 
the domain of asymptotically high energies. The experiments at RHIC provide a powerful argument that the domain 
of interest in heavy-ion collisions is a strongly coupled domain. It will be very interesting to see to what extent this 
remains true at LHC energies. 

Also notable in this regard are lattice calculations of the spectral function. The striking feature of the lattice results 
is the sign oscillation of the spectral function in the w ~ 2'kT region. This oscillation is absent in weakly coupled 
calculations [U,!!^, but appears to be a generic feature in the spectral functions obtained in holographic models, as 



our results in Figs. 1(b) and [5] illustrate. This qualitative difference appears to be another manifestation of the now 
familiar fact that real-time response (in particular, hydrodynamics) is more sensitive to the coupling strength than, 
e.g., equation of state. 
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Appendix A: Background of the holographic model 

1. Einstein's equations and one-point functions 

This appendix summarizes, for completeness, the relevant results of Ref. [2^. The background geometry corre- 
sponding to homogeneous boundary conditions on g^i, and (j) can be found by solving Einstein's equations for the 
metric in Eq. which take the form: 

B = -\^^ (Al) 
/ = (4 + b) /, (A2) 

-6/ -I- / (24 - + 2e2-B V{4>) = 0, (A3) 

0/ + - /(4 + B)) - dV{cl,)/dcf> = 0, (A4) 

where a dot denotes a logz derivative, e.g., (j) — zd(j)/dz. The holographic correspondence provides the boundary 
conditions at the UV boundary z = e given by Eqs. ((23)) . Minkowski metric at the boundary requires 



fie) = 1 . (A5) 

Equation (jA2p can be integrated once to give 

/ = -wz*e^. (A6) 



^ Holographic models which reproduce the effect of the logarithmic running of as, have been discussed in, e.g., Refs.[45l447|. 
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The integration constant w must be positive if the metric is to possess a horizon /{zh) = at some value of zh- 
Since w determines the position of the horizon, it is related to temperature, and we find that it is proportional to the 
enthalpy, Eq. (jA12j) . 

A boundary condition on B is not needed because the value of B is determined by Eq. (IA3p . which is algebraic 
in B. The role of the second boundary condition for Eq. (jA4p is played by the requirement that is finite at the 
horizon, z = zh, which is a regular singular point of the second order differential equation (|A4p . 

Near the z = e boundary, 0— >-0, i?— )-0, 5— >-0 and / — 0. Equation (jA4l) for (j) can be linearized and the 
asymptotic behavior of near the boundary can be determined easily: 

(/.(z) ^{c- de'^+~'^- ) z-^- (1 + . . .) + (1 + . . .), (A7) 



where the curvature of the potential V"{0) = m| determines the indices A± — 2±y/4 + -m?. The coefficient of the first 
term is related to c by Eq. The coefficient d of the second linearly independent solution should be determined 

by the finiteness condition at the horizon and is a function of w (i.e., temperature) and c. 

By calculating the derivative of the 5D action with respect to c and matching it, by holographic correspondence, 
to the expectation value (O), one finds (see also [48|) 



dc 2k2z3-A- 



where ". . . " denote UV divergent but temperature independent terms. One can thus see that the integration constant 
d is related to the UV finite (and temperature dependent) part of (O). 

Furthermore, from the expression for 0(z) near the boundary Eq. (jA7l) and Eq. (jAip . the function B{z) can be 
calculated for small z: 

B{z) = zB'{z) = -^-Al (c - de^^-^-f z^A- „ Ia^ A+ (c - ds^+-^-) dz^^+^- .... (A9) 
6 3 

By considering homogeneous variations of the boundary condition on the metric, the one-point functions of the 
stress-energy tensor can be calculated using Eq. ^ as in Ref. ^29|] . 



6e-s(z) 



(r") = ^-(r™). (Aio) 



2k2z4 

The thermal energy and pressure, 

6 - (TO") - (TOO)t=o ^ [(T"°)] ^, p = (T") - (T")t=o ^ [(T")] ^, (All) 
are finite at e = and equal to zero at T = 0. The enthalpy, w = e + p is related to the constant w: 

(A12) 

After solving Eq. (jASp for i? at z = £ with given by Eq. (jA7p , we find that the energy and pressure can be expressed 

as 

^=f-^¥^-(^+-^-)' P = ^-e. (A13) 
4 8k'' 

Therefore the temperature dependence of the expectation value of the trace anomaly is given by, 

= 3p - 6 ^ (A+ - A_) = -A_ c [(O)]^. (AM) 

in accordance with the anomaly equation (jl9p . 



2. High temperature limit 



The equations of motion can be simplified and analytically solved if one considers the high temperature limit, or 
the limit where T ^ c^l^- . Since the enthalpy is related to the temperature, 

w = 4(^r)*(l + 0(c2T-2A-))^ (A15) 
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this limit can also be expressed as w ^ c^l ^ - . One can begin by observing that at large w the function / varies very 
rapidly according to Eq. (|A6p . This means one can neglect variation of the function B between the boundary z — e 
and the horizon z = zh^ since zh becomes small (as ui~^/^). Since on the boundary _B = (up to terms of order 
£^^- , negligible here, according to Eq. (jASp ). we find from Eq. (|A6|) 



/(z) = l-u;z^/4. (A16) 

Another consequence is that which is small at z = e, remains small up to zh {4> ^ cz^" ^ cjT^^ <C 1), and the 
linearized approximation to Eq. ()A4|) is valid not only near the boundary, but all the way to the horizon. With B = 
and / from Eq. (|A16|) we obtain 

l-l..')*"-(f + l^f),'-!^*.0^ (A17) 

Equation (jA17|) can be solved analytically 

(t,{z) = cz^- 2i^i (A_/4, A_/4, A_/2, «;z74) 
+ dz^+ 2i^i (A+/4, A+/4, A+/2, u;zV4) , 

where the coefficients follow the notations of Eq. (jA7p (up to terms 0{e^+^^-), here negligible). Both linearly 
independent solutions are logarithmically divergent at the horizon z — zh, where wzjj/A = 1. The condition |0(z^f)| < 
oo requires us to select the linear combination in which these divergences cancel. This fixes d in terms of c: 

d=-ci(;(^+-^-)/4i?(A_), (A19) 

where the function D(A_) = 1/D(A+) is given by 

Z5(A_) = Z^cot(7rA_/4)4^;^. (A20) 
V y 2- A_ ^ ' ' r(A_/4)4 ^ ' 

This solution for (p{z) can be used in Eq. (|Aip to calculate B'{z), 

i3'(z) = -i(cA_z^-+i2Fi(A_/4, l + A_/4, A_/2,u;zV4) ^ 
6 (A21) 

+ dA+z'^++\Fi (A+/4, 1 + A+/4, A+/2, wz^/A)f 

From this expression, the function B{z) can be calculated to order c^, which is the leading order in the high- 
temperature limit. By iteratively solving the equations of motion, the higher order corrections to B as well as to / 
and (j) can be determined if necessary. 
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